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Abstract 

We investigate which experiments are better suited to test the robust prediction that cos- 
mic strings do not produce secondary Doppler peaks. We propose a statistic for detecting 

oscillations in the C' spectrum, and study its statistical relevance given the truth of an infla- 
tionary competitor to cosmic strings. The analysis is performed for single-dish experiments 
and interferometers, subject to a variety of noise levels and scanning features. A high resolu- 
tion of 0.2 degrees is found to be required for single-dish experiments with realistic levels of 
noise. Interferometers appear to be more suitable for detecting this signal. 
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In two recent Letters |^ it was shown how the existence or absence of secondary Doppler 
peaks in the cosmic microwave background (CMB) power spectrum C' could rule out or confirm a 
large class of defect theories, including cosmic strings. The argument used is attractive in that it 
does not depend on details of defects or inflation, but only makes use of well-established contrasting 
properties peculiar to the two types of theory. This is particularly welcome when calculations in 
defect scenarios are so difficult and unreliable. A simple but robust test based on an issue about 
to be decided by experiment seems a soundly cautious approach to defect theories. 

The question therefore arises: which CMB experiments can resolve the secondary Doppler 
peaks? This is a timely issue when so many proposals for ground based and satellite borne 
CMB experiments are being made [H . Experimental features have so far been motivated by their 
implications on inflationary parameter fixing from Doppler peaks' position, height, and shape Q |^ . 
Secondary peak detection is a far less demanding task, and can be used to quantify experimental 
spectral resolution at the most basic level. In we analyze this problem in a general setting, but 
here we consider only its bearings on cosmic string theories. 

Apart from the absence of secondary Doppler peaks in cosmic string theories, the only reliable 
feature known is that their primary peak is located at I = 400 — 600. If the main peak is measured 
outside this range one can rule out cosmic strings, but we shall assume here that this is not the case. 
In order to identify experiments targeting the string's lack of secondary oscillations we investigate 
how we could falsify cosmic strings given the truth of a competitor inflationary scenario with a 
primary peak located in I = 400 — 600 and exhibiting secondary oscillations. For definiteness we 
have taken CDM with Q, = 0.3, h = 0.6, ^bh? = 0.02, and a fiat primordial spectrum. We shall 
call this theory stCDM, the CDM competitor of cosmic strings, and we plot its C' spectrum in 
Fig |l|. The main peak height and shape will be assumed to be the same for cosmic strings and 
stCDM, and the low I section of the spectrum will be ignored. In this way we assume maximal 
confusion in whatever is uncertain, or alien to the signal to be experimentally tested. 

The idea is to apply to stCDM a statistic sensitive only to the existence or absence of secondary 
oscillations. To set it up, we first compute the average power Ci in bins i = 1,2,3 denoted by 
horizontal bars in Fig [l[ We then infer the spectrum convexity with C — (Ci + C3)/2 — C2. 
A positive convexity reflects unambiguously the flrst dip of the spectrum and therefore the flrst 
secondary oscillation. If the overall error in C is cr(C) then one can claim that C is positive (and 
therefore that there are secondary oscillations) with a number of sigmas equal to 

S is then to be seen as the stCDM secondary peak detection function, or equivalently, the cosmic 
string rejection function. 

In this Letter we set up a large parameter space of experiments, on which we compute the 
contours of S. We consider two types of experiments: single-dish experiments (recovering some 
of the results in ||^) and interferometers. For single-dish experiments we allow the beam size, sky 
coverage, and detector noise to vary. For interferometers we take as free parameters the primary 
beam, number of fields, and detector noise. We consider errors associated with cosmic/sample 
variance, spectral resolution limitations due to finite sky coverage, and instrumental noise. Fore- 
ground subtraction errors are included (not naively, as we prove in in the form of only an 
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extra instrumental noise term. 

We outline a method for computing errors in estimates explained in more detail and gen- 
erality in For simplicity we consider the small field limit, and assume at first no instrumental 
noise or foreground subtraction uncertainties. We stereographically project the sky onto a plane, 
and expand in Fourier modes, using the symmetric notation (factors of 27r evenly distributed). We 
denote by a(k) the modes provided by an all-sky observation with infinite resolution, and a''(k) 
the modes as seen through an observation window W{'x.) and convolved with a beam -B(x). For a 
single-dish experiment we shall assume that the window is a square of side L treated with a cosine 
bell 1^] (to bar edge effects), and that the beam is a Gaussian with FWHM 0^. For interferometers 
the window (better known as the primary beam) is a Gaussian with FWHM 6^ , and the beam is 
essentially unity j|] . Using the convolution theorem twice we have that 

a^(k) = J ^a{k')B{k')Wik k') (2) 

where W^(k) is the window Fourier transform. The all-sky modes form a diagonal covariance 
matrix (a(k)a*(k')) = C(A:)5(k— k'), where the brackets denote ensemble averages. In calculations 
concerning small patches of the sky C{k) can be obtained by interpolating the with k = I 
On the other hand the sampled modes a*(k) form the covariance matrix |^ 

(a'*(k)a^*(k')) = 

/ ^C{k")B^(k")W{k - k")W^*(k' - k") (3) 

which encodes all the finite sampling hurdles for recovering the power spectrum C', now to be 
examined. 

Firstly, the sampled modes power spectrum C''(k) = (a*(k)a**(k)) becomes the convolution 
of the raw power spectrum with the window power spectrum ||, ^ ^ . This has the effect of 
leaving a low k white noise tail in C"'(k) up to fc « 1/L, and thereafter averaging out oscillations 
in the raw spectrum on a scale Al ^ 1/L. If the field has edges there will also be spurious 
oscillations of period 1/L superposed on the spectrum. Field edges can be treated as in §. In the 
case of a square field the window should be multiplied by a cosine bell. Whenever the sampled 
spectrum is highly distorted, a deconvolution recipe is then required. In the presence of noise 
and cosmic/sample variance this induces a large deterioration of the detection function. We have 
however checked [|| that, providing edge effects are treated, one has in the stCDM Doppler peak 
region 

/lik' 
^\W{k')f^aC{k)B'{k) (4) 

for fields with i > 4 degrees, or 0^, > 2 degrees (we illustrate this point in Fig. [l]). Therefore, 
as long as the field is not too small, deconvolution is trivial in the relevant spectrum sections of 
stCDM, and it does not add any extra errors. Note that a = 0*/(27r)^, where il^ is the field area 
(O^ = I/^ for a square field, fi^ = tt^^ / (8 log 2) for a primary beam) . 

Secondly, finite fields have the effect of correlating neighbouring a" (k) modes within a corre- 
lation radius ^ w 1/L (assuming edge effects have been treated) or ^ w l/^ui- This translates 
into a correlation length ^(fc) between C{k) estimates. In M we show how spectral resolution 
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Figure 1: The angular power spectrum of stCDM (line) and one possibility for cosmic strings 
(points). The dashed curve is the stCDM power spectrum C*(fc) as sampled by an interferom- 
eter with 9w = 2° (divided by a defined in Eqn. (||)). We have plotted the bins used in C as 
horizontal bars, their level indicating the average spectrum in bin. The vertical errorbars are the 
cosmic/sample variance. 

in the Doppler peak region is typically imposed, not by the fact that a given C'{k) estimate re- 
ceives contributions from neighbouring fc, but because we can only make uncorrelated estimates 
of the power spectrum with a separation ^(fc). This effect is reminiscent of cosmic covariance 
in non-Gaussian theories p^ . Correlations also determine the cosmic/sample variance. Using 
cov{x^ ,y^) — 2cov^{x,y), it can be proved that any power spectrum estimate Cn (using a 2D 
region fl of the Fourier domain, with area Ak, in which C"*(k) does not vary much) is affected by 
the sample variance 



Nq acts as the effective number of independent modes in the region, and it can be used to define 
an average density of independent modes po- 

We have found it convenient to replace the k space by a square mesh, to be called uncorrelated- 
mesh, with a spacing locally given by feg ~ 1/VPo- This mesh, on average, contains all the 
non-redundant information, given cosmic/sample variance and the correlations imposed by finite 
sampling. We have checked that the uncorrelated-mesh is nearly a square lattice with ko « 27r/L 
for a square field, and fcp « 2\/47r log 2/0^ for a Gaussian field. Although it is easy to improve on 
this approximation, it is normally a good enough recipe. Using only mesh points (denoted by k^) 
the sampled power spectrum C^{k) can be estimated with 




(5) 




(6) 
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where Nk is the number of modes in the mesh which satisfy |ki| = k. The residual correlations 
between these estimates fall below 5%, but only a finite number of k can be estimated. Their 
average separation Ak ~ -y/fc^ + /cq/tt — k, for k > ko, defines the maximal spectral resolution. 
More estimates could be inserted in between these, but they would necessarily be correlated 
estimates. Only for k > k'^/{2TT) can individual C' be estimated (A/c w 1). The cosmic variance in 
these estimates is approximately cr^(Cs(fc)) w 2C'^^{k)/Nk- For k > kQ/{2n) this means a^{C^) ~ 
(C'V0(47r/L2), as naively expected For k < kl/{2iT) the naive expectation breaks down. 

We now study the effects of noise, differentiating between single-dish experiments and interfer- 
ometers. Let Q.pix be the pixel area, and cr^^^ be the noise per pixel [Q. Fixing the detector sensi- 
tivity s and total time of observation ttot fixes the quantity w"^ = Airs'^/ttot = fp^^pia; (47r/L^), 
which we therefore use to parameterize the noise level. The noise introduces an extra diagonal term 
with value aw^^L'^/{AiT) into the mesh modes covariance matrix. Hence a centred uncorrelated- 
mesh estimator of the power spectrum is now 



\ki\=k 



and its variance is 



g'(c(fc)) _ 2 / w-^L^ y 

^ Nk\ ^ 4TTB^C{k) J 



(8) 



For interferometers p4| ] noise is added directly in Fourier space. The noise in a given mesh cell 
is given by aj^ = s^Vt'''^ /{tvisNms)i where Nyis is the number of visibilities in the cell, s is the 
sensitivity of the detectors, and tyig is the time spent observing each visibility. The coverage 
density pc — N^stvis^^ /tf (where i/ is the time spent on a given field) is assumed to be uniform 
in a ring of the itw-plane containing the stCDM relevant bins. This assumes a dish geometry 
like the one proposed in (iH). If one decides to observe n/ well-separated fields, then each mesh- 
point acquires an extra index i = I, . . . ,n/, and points with different indices are uncorrelated. 
For fixed detector sensitivity and total observation time one should now keep constant = 
{2itY / [pcttot) = (27r)^(T^/(ri^^n/), and so this is the noise parameterization we choose for 
interferometer estimates. A centred estimator is now 

cwH^Ei'-iMM-^ (9) 



and its variance is 



|k,|=fc 



(T^{C{k)) _ 2 ( w-^Q'^Uf 



(10) 



C2(fc) Nk V C{k) 

From these results one can compute E in the large experiment parameter space proposed (which 
always assumes L > 4°, or 9^ > 2°). Two types of results ensue. Firstly, one can provide guidance 
on experimental design given a constraint, such as finite funding. This constraint is mathematically 
translated into hypersurfaces of constant w^^ . Secondly, we may provide the value of the detection 
as a function of w^^, assuming ideal design. This will set lower bounds on w^^ for any meaningful 
detection, telling us thereafter how fast the detection improves with a given w^^ improvement. 
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Figure 2: Low noise [w ^ ~ (25^ii')^(°)^, top) and high noise {w ^ — {60^K)'^{°)^ , bottom) 
contours of the 'S{L,9i,;w^^) function for stCDM. We have also plotted the ideal scanning lines 
i.(0fa). 



In Fig. 1^ we show a low and a high noise section of 9i,; w^^) for single-dish experiments 
p6[ . Most noticeable are the high resolutions required for a significant detection {9b < 0.3° and 
9i, < 0.25°, respectively). These are due to the fact that we are testing features at a rather large 
I, and the noise term goes up exponentially with I as we approach the resolution limit. It is 
also obvious that all-sky scanning is not ideal under realistic levels of noise. For fixed resolution 
and noise, increasing the coverage area will at first increase the detection, but beyond a certain 
coverage Li, the detection will initially saturate, then start to decrease. This is because noise 
separation relies on allowing the same coherent signal to compete with the incoherent noise. Only 
after the signal has dominated the noise does it make sense to increase the coverage area, so as 
to reduce the sample variance. If the noise is very high, then all Uot should possibly be used in 
a small patch of the sky (larger than 4°). The ideal scanning lines Li{9b) are plotted in Fig. 
As the resolution increases so do Li and the achieved Ti{Li,Ob\w~^). Initially they increase very 
fast; then, for 9^ < 0.1°, not by much. For a low noise experiment {w^^ — (25/ii4r)^(°)^) the 
detection increases from S = 1 to E w 36 as the resolution is improved from 9i, — 0.3° to 9^ — 0.1 
(with Li 4° and Li « 65°, respectively). From then on S improves by only a few sigmas. 
Even with infinite resolution the ideal coverage area is L « 80°. For a high noise experiment 
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Figure 3: ^{0w — 2°,n/; w^^) contours and density map, and ideal scanning line nfi{w~^). 

{w^^ ~ {60 jj.K)^{°)^) the resolution is even more crucial. One needs 6b = 0.25° to obtain a 3 
sigma detection (with L = 4°), and a decent E = 10 detection can be achieved only with 6'f, « 0.15° 
(and Li — 18°). It can be checked that all-sky coverage is only useful if the noise is lower than 
about = (11^X)2(°)2. 

In Fig. ^ we show l](0iu = 2°, n/; w^^) for an interferometer. Ideal scanning now always means 
9m = 2°, and the ideal coverage area is fixed by a curve nfi{w~^). The high noise region is 
> {150fiK)'^iad^^. There one should look into only one or two 2° fields, in order to obtain 
a detection between 1 and 2 sigma. For w^^ < (150/iif )2rad^^ we enter the signal dominated 
region. Following the ideal scanning line with decreasing w^^, Ufi and S start increasing, first 
slowly, then very fast. For w^'^ = {lOQ^KYrad"^ one may obtain a 3 sigma detection using 8 
independent fields with 9^ — 2°. For w^^ = (20/j,i^)2rad~^ it is worth looking into about 40 
of these fields, obtaining thus an 8 sigma detection. We have estimated CAT noise level to be 
w^^ = (300/iiir)2rad^®. This is a mere prototype, and a 10-fold improvement should be easily 
attained. 

These results stress the contradictions of an all-purpose experiment. If the low-^ plateau of 
the spectrum is the theoretical target then one needs all-sky coverage, and satellite single-dish 
experiments are to be favoured. Even if one wishes to target the main sCDM features, encoded 
mostly in the first peak's place and height, then this is still true ||l^. Our work shows how such 
a design relies heavily on the assumption that the signal is in the vicinity of sCDM. If instead 
one is to test the high-Z opposition between low il CDM and cosmic strings, then we have seen 
that single-dish experiments are required to have rather high resolutions. Interferometers appear 
to be less constrained, providing 2-3 sigma detections under very unassuming conditions, with 
rapid improvements following further experimental condition improvement. Furthermore, in this 
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context, all-sky scanning is not only unnecessary, but in fact undesirable. The best scanning 
is normally achieved with deep small patches. These two features contradict sharply the ideal 
experimental design motivated by the standard theoretical gospel. We believe that a variety of 
contrasting experimental techniques may equally well find their niche of important theoretical 
implications. 

We should mention, in closing, that if one is to combine the high-^ cosmic string signal with the 
requirement that the low-/ section of the spectrum is to be measured, then the logic is naturally 
changed. See |0 for this alternative perspective. 
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